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We use the identification of the edge mode of the filling fraction v = 1 quantum Hall phase with 
a 1+1 dimensional chiral Dirac fermion to construct an analogue model for a chiral fermion in a 
space-time geometry possessing an event horizon. By solving the model in the lowest Landau level, 
we show that the event horizon emits particles and holes with a thermal spectrum. Each emitted 

O 

quasiparticle is correlated with an opposite-energy partner on the other side of the event horizon. 

U 

Once we trace out these "unobservable" partners, we are left with a thermal density matrix. 
> PACS numbers: 04.20.-q, 04.70.-s, 73.43.-f 
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I. INTRODUCTION 



There are several apparently different explanations for the origin of black hole radiation. 
In his original account pQ Hawking kept track of what one means by a "particle" as a 
wavefunction propagates in the background geometry. A field theory derivation using the 
trace anomaly in the energy momentum tensor was given given by Christensen and Fulling 
[2], and more recently Robinson and Wilczek [3] and others [U [5] have applied the two- 
dimensional gravitational anomaly in the region near the horizon. Yet another route obtains 
the Hawking radiation from quantum tunnelling across the horizon [Bj. (For a review of the 
tunnelling approach see [7].) 

Given these alternative derivations, it is reasonable to ask just what is required for an 
event horizon to emit thermal radiation. Is gravity really necessary? This question has led to 
the study of analogues of black holes and event horizons in other areas of wave propagation. 
The first such analogue was the acoustic black hole proposed by Unruh, who discovered that 
the wave equation for sound in a background fluid flow was equivalent to the wave equation 
for a scalar field in a curved space-time [8]. The subject has now developed extensively, with 
gravity and Hawking radiation analogues being proposed and constructed in quantum-fluids, 
optics, and solid-state devices. For review with an extensive list of references see [9]. 

The present paper proposes a conceptually simple, and possibly experimentally realizable, 
condensed matter model of quantum mode propagation in which an event horizon emits 
thermal radiation. The analogue space-time is flat, but consists of two causally disconnected 
halves. It is therefore a member of the general class of condensed-matter event horizons 
discussed by Volovik in [TO]. Our model exploits the intepretation of the edge-modes of a 
filling fraction v — 1 quantum Hall system as a massless chiral Dirac fermion whose local 
"speed of light" is determined by the potential that confines the Hall fluid, and is therefore 
subject to external control. 

In the next section we describe the model in the language of first- quantized tunnelling. 
In the third section we adopt a second-quantized formalism so as to obtain a Bogoliubov 
transformation between two natural bases for the system. This allows us to display the 
physical "vacuum" as a coherent superposition of particle-hole pairs that are entangled 
across the horizon. Just as the Minkowski pure-state vacuum is a thermal mixed state when 
seen by a Rindler co-ordinate observer [TTjlT2"]. our pure-state vacuum appears thermal when 
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we trace out the "unobservable" over-the-horizon member of each pair. 

II. LOWEST LANDAU LEVEL 

We model our black hole as a two-dimensional electron gas (2DEG) in the v — 1 quantum 
Hall phase. We arrange for the boundary of the 2DEG to lie along the y axis, with the region 
x < occupied by the gas, and the region x > empty. Now assume that we have engineered 
the "confining" potential V to be of the form 



and have chosen the the direction of the perpendicular magnetic field B so that the classical 
guiding-centre drift velocity is 



The electrons move along equipotentials V(x, y) = E, which, for this potential, are rectan- 
gular hyperbola? that have the x and y axes as asymptotes. In particular, the electrons at 
the edge of our 2DEG move vertically along the y axis at velocity 



This velocity is our analogue of the local speed of light. The edge modes in the regions y > 
and y < move in opposite directions, and so these two regions are causally disconnected. 
They are separated by an event horizon at y — 0. 

The price we pay for the event horizon is that the electrons in the occupied region with 
y < (the interior of the black hole) are in a state of population inversion. In the absence of 
the magnetic field the electrons would rapidly fall into one of the the lower energy quadrants. 
Because of the strong field, however, and in the absence of inelastic or tunnelling processes, 
they are constrained to stay on their hyperbolic classical orbits. The inherent instability 
of the "vacuum" in the black hole interior corresponds to the observation of Parikh and 
Wilczek [6j that a black hole must be thought of as highly excited quantum state. 

Except for the case E = 0, each of the classical equipotentials Xxy = E consists of two 
disconnected branches and intitally all the particles lie on only one of these branches. The 



V(x,y) = Xxy, 
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FIG. 1: The 2DEG black-hole analogue. The shaded region is the 2DEG. The lines indicate the 
semiclassical electron orbits (dashed when mostly unoccupied). The low energy excitations near 
the boundary at x = constitute the quantum system in which we will find Hawking radiation. 
This radiation is illustrated by three correlated pairs of electrons and holes moving in opposite 
directions inside and outside the black hole. 

branches for small E approach each other near the origin. There is therefore a non-zero 
amplitude for a particle to tunnel from one branch to the other of the same energy. This 
tunnelling leads to electrons and holes being emitted from the event horizon. 

To calculate the tunneling amplitude, we will assume that the magnetic field is large 
enough that we can ignore all Landau levels except the lowest. The lowest Landau level 
(LLL) approximation is very natural as it is this situation that the excitations near the edge 
of a quantum Hall droplet can be identified with those of a 1+1 dimensional chiral fermion 
with Hamiltonian 

/oo 
v e d g c{y)4> ] {-idy)i)dy. (4) 
-oo 

In this picture the 2DEG itself is the filled Dirac sea. 
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We chose the symmetric gauge in which the LLL wave-functions are of the form 

ip{x,y) = exp j-^-e^l^fj^tz), (5) 

where z = x + iy. All quantum information resides in the holomorphic factor ip(z), and we 
will refer to this factor as the LLL "wave-function." We therefore regard the LLL Hilbert 
space as a Bargmann-Fock space of finite-norm holomorphic functions with inner product 

(<p,X) = J d 2 ze- eB \ z \ 2/2 ^{z)x{z), d 2 z= j.dzAdz = dxAdy. (6) 

Bear in mind however that the LLL wavefunction should be multiplied by exp { — |e5|2;| 2 } 
before plotting probability densities or computing currents. 

The action of z on the LLL wavefunction is by simple multiplication, but multiplication 
by z takes us out of the space of holomorphic functions. The LLL operator corresponding to 
z becomes instead z\ where the adjoint is taken with respect to the Bargmann-Fock inner 
product. This identification makes 



For our potential 
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Z ^ z] = jBdz- (7) 



Xxy=^(z 2 -z 2 ) (8) 



the first-quantized eigenvalue problem 

Htf) = eip (9) 

therefore becomes 

.J*j? ~ t) !{:) - (10) 

Only the potential appears in the this equation as the LLL wavefunctions are annihilated 
by the electron kinetic energy operator. A rescaling gives us a standard form of Weber's 
equation (See [13J §16.5, or [TJ] chapter 19.): 

d 2 r 



d( 2 4 



/(C) = o/(C), (ii) 



with 

'eB 



a = -te[ — )f(C), C = VeBz = —. (12) 
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FIG. 2: Left figure: A density plot of the the absolute value of the even function 
exp{ — \z\ 2 /A}yi(x, y) for the case e = — 10. Right figure: A density plot of the absolute value 
of the odd function exp{— \z\ 2 /4}y2(x, y) for e = —2. 

For simplicity we now set A = eB = 1. We can always restore the general parameters by 
scaling the units of length and energy. 
If <f(e, z) is a solution of 

(|^-£)^C) = -^(e,C) (13) 

then so are <p(e, —z), <p(—e, iz) and ip(— e, — iz). At most two of these solutions can be linearly 
independent. 

A fundamental pair of independent solutions is 



iFi U 



e 1 z 2 
i 

3 e 3 z 2 



V2(e,z) = ze -'/\Fi{--i-,-,^y (14) 

Here i-Fi(a, b, z) is the confluent hypergeometric function. These functions are even and odd, 
respectively, under z —z. After multiplication by exp{ — 1^| 2 /4} the resulting wavefunc- 
tions are localized on the semiclassical orbits, which form the two disconnected branches of 
the rectangular hyperbola xy = e (see figure [2]). These solutions to the LLL potential have 
been studied in connection with Riemann hypothesis |15j . 



More useful to us is the solution of ( 13 ) given by the parabolic cylinder function 



2 -(fc/2+V4)_l 

— \/2 sin 



Hi 



7T 



r 



cos 

1 

4 



7T 



ze 
2~ 



ze 
1 



ze 
2 



I _I+2 

2 e 2 



r 



3 ze 
i + 2 



(15) 



,2 ^0 

Here D n (z) is Whittaker and Watson's notation for their parabolic cylinder function 
and U n (z) is the now more common notation used by Abramowitz and Stegun |14j . The 
essential properties of U n {z) are that it is an entire function, and that it decays rapidly as 
x — > +00 for any real or complex n. 

The solution /7_j e (z) describes particles moving in from the left (the occupied region) in 
the lower left quadrant if e > and the upper left quadrant if e < 0. They mostly remain 
in that quadrant, but there is some probability of tunnelling to the other branch of the 
hyperbola (see figure [3]). If e > the result is that a tunnelled positive energy particle 
is emitted by the black hole, leaving a negative energy hole (i.e. the absence of positive 
energy particle) inside the event horizon. If e < then a positive energy hole (the absence 
of a negative energy particle) is emitted by the black hole leaving a negative energy particle 
inside the event horizon. 

Along with the solution U_ ie (z) we have the solutions LLj e (— z) and U ie {iz) and U ie (—iz). 
We will find use for all of these solutions, as they describe motion with different boundary 
conditions (see figures [i] and [5]) . 

To discover the relative amplitudes of the direct and tunnelled waves we can use the 
asymptotic expansion 



Near the y axis this reduces to 



Z ) rv, e ~\z\ 2 /4-z 2 /4 z ie-l/2 



O 



axg(z)\ < 3/4. 



(16) 



if)(x,y) ~ (gauge phase)e 



-x 2 /2 



expjieln \y\ — sgn (y)e7r/2}. 



(17) 



The ratio of tunneled to direct amplitude is therefore exactly exp{— ne}. 
We can conform this result by using the identity 

r(| + ze) 



UJz) 



271 



[e~ en/2 e- ln/i U ie {tz) + e e7T/2 e +i7T/4 U it {-tz)} 



(18) 
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FIG. 3: A countour and a density plot of exp{ — \z\ 2 /4}U-i e (z) for the case e = — .5. Particles enter 
from the left and the beam divides between down-going and weaker tunnelled edge-mode wave and 
the up-going and stronger direct edge-mode wave. 

together with the fact that LLj e (z) tends rapidly to zero in the right half-plane. Thus, if R 
is positive 



U-u{iR) 




[e- €7T/2 e- in/4 U ie {-R) + e +OT/2 e +w/4 f4(i?)] 



r(| + *Q c 




/2 e- i7v/ *U ie (-R) 



(19) 



and 



U-u{-iR) 



r(i+ig) 

r(i + xe) 



[e tn/2 e- ln/i U lt (R) + e tn/2 e +i ^U ie (-R)] 




(20) 




The direct and tunneling amplitudes therefore have magnitude 





(21) 



Note that \d(e)\ 2 + \t(e)\ 2 = 1 because r(z)r(l 



z ) = 7rcosec (nz ) gives us 




(22) 
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The occupation probability of an outgoing particle or hole state with energy e is therefore 

Pie) = — f 7 . (23) 

W 1 + exp{27re} V ; 

The chiral edge states emerging from the event horizon are therefore thermal with 

T =h < 24 > 

or 

k ° T = < 25 > 

once we restore parameters and units. Comparison of this with the usual Hawking radiation 
formula 

^B^Hawking = — (26) 

indicates that the "surface gravity" k of our analogue black hole is the edge-velocity accel- 
eration 

(27) 



eB dy 



horizon 



III. SECOND QUANTIZATION, MODE EXPANSIONS, AND A BOGOLIUBOV 
TRANFORMATION 

The space of LLL functions ^ does not contain the delta function. Its place is taken by 
a reproducing kernel 

r t def 1 [ 1 . ,n 1 I io / 1 _ I 

{xt^yx^y^} = —exp<--\z 1 \ - -\z 2 \ /4 + -z x z 2 > . (28) 

If ip(x, y) is of the form ^ then 

J d 2 z 1 ?p(x 1 ,y 1 ){xi,y 1 \x 2 ,y 2 } = 4>(x 2 ,y 2 ). (29) 

In particular {xo, yo\xi, yi}, considered as a function of (xi, yi), is of this form, so the kernel 
reproduces itself: 

J d 2 z 1 {x ,y \x 1 ,y 1 }{x 1 ,y 1 \x 2 ,y 2 } = {x ,y \x 2 ,y 2 }. (30) 
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When we expand out the second quantized LLL field operator in terms of the discrete set 



of normalized eigenmodes z n / \j2n2 n n\ for the potential 



we find 



V(x,y) = l -(x 2 + y 2 ) = ^zz^z^ (31) 



oo 1 

~ V2n2 n n\ 



n=0 

Here the operators a n obey 

{a n , & m } — 5 nm , (33) 

and the usual canonical anticommutation relation {ip^(x),ip(x')} = 5(x — x') for the field is 
replaced by 

{^{x 1 ,y 1 ),^{x 2 ,V2)} = {x 1 ,y l \x 2 ,y 2 }- (34) 
If we retain only the holomorphic factors, then we have 

We can also expand in a continuous set of continuous set of eigenfunctions. For example, 
we can make use of energy E eigenfunctions for the potential V(x, y) = x. These are 

M*) = ^71 ex P [ Ez - -/ ~\&\- ( 36 ) 

They have been normalized so that 

{<p E ,<p E ,)=2n5{E-E?). (37) 
The holomorphic field operator is then 

/OO 1£J 
— a E ip E (z) (38) 

with 

{a E ,a j E ,} = 2n5{E - E'). (39) 



We easily confirm that (38) still satisfies (35). 

Similarly, we can expand the LLL field operator in terms of complete sets of parabolic 
cylinder functions. There are two distinct ways of doing this. Begin by defining 

<p^ + \z) = ±r(l/2 -ie)U-*(z), (40) 



V^'K*) = ^74 r (V2 - ie)EU(-*). (41) 
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a) 



b) 



FIG. 4: The "in" wave functions: a) (fi in ' + \z) for e<0,b) for e > 0, c) ^ in ' \z) for e > 0, 



d) ipi m ' \z) for e < 0. In each case the incoming wave divides between two outgoing waves. 



The label "in" designates that these wave-functions describe states that have a simple de- 
scription prior to their opportunity for tunnelling (see figure [1]) . These "in" functions have 
been normalized so that 



(cp^ a \^' a,) ) = 2n5(e - e')5 aa ', a = ±, 



(42) 



and they obey the LLL completeness relation 

E 



a=± 



fc fin, a)/ \ (in,«)/ \ 1 11 



2tt 



2,7i 6XP 2 ZlZ > 



(43) 



(Both normalization and completeness are easily established from the integral expression in 



the last line of (15).) Then we can set 

de 



2tt 



(44) 



The "in" vacuum is the appropriate many-body state for our initial conditions. It is charac- 
terized physically by the condition that no particle is approaching the 2DEG from the empty 
single-particle states to the right, and that all the single-particle states incoming from the 
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a) 



b) 



c) 



FIG. 5: The "out" wave functions: a) ^ ou *- ert ) (z) for e > 0, bj ^ out ' ext) (z) for e < 0. c) ^° ut,int) (z) 
for e > 0, d) (pi° ut ' mt \z) for e < 0. In each case two weaker incoming waves assemble the outgoing 
wave. 



left are occupied. It is characterized mathematically by the conditions 

a e |0,in) = = 6 £ |0,in), Ve. 



(45) 



The second set of functions is 



^(out.cxt)^ 



71-1/4 



F{l/2 + ie)U ie {iz) 



¥> ( r tM) (z) = ^(1/2 + ie)U i£ (-iz). 



They are also orthogonal 



(out, a) (out, a') 



) = 2n 5(e - e')5 aa ' . a = int, ext, 



and obey the LLL completeness relation 

de 



2tt 



<P. 



Z2) = 2^ eXP ]2 ZlZ2 



(46) 
(47) 

(48) 

(49) 



The labels "ext" and "int" indicate that the functions live mostly in the exterior (y > 0) 
and interior (y < 0) of the black hole. They decay rapidly in the other region (see figure [5]). 
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In terms of these new functions we have 



— a<? ut ' a V° nt,a) (z) 
00 27T 



(50) 



The "out" operators a^ 1 '^ and (a^ ut ' a ^y create and annihilate particles that are simply 
described as excitations over the asymptotic naive vacuum in which every state in the region 
x < is filled and every state in x > is empty. For e > the operator o;( out ' cxt ) annihilates 
a positive energy particle in the asymptotic region y ^> outside the black hole. For e < it 
annihilates a particle in the 2DEG and so creates a positive energy hole in the same region. 
For the a( out,mt ) that act on states within the black hole the roles of hole creation and particle 
annihilation are reversed as the 2DEG consists of particles with positive energy. To stress 
the causally disconnected character of the interior and exterior regions, we will write "out" 
vacuum as 

|0,out) = |0,out,ext) <S> |0,out,int) (51) 

with 



2<°uW)|0,out,ext) = 0, e > 

^(out,cxt)) t | 0)Out?ext } = 0> e < Q) 



(52) 



and 



a(° ut ' int) |0,out,int) = 0, e < 
^(out,i„t))t| 0?out)irit ) = 0, e > 0. 

Comparing the two expressions for ip(z) gives us the Bogoliubov transformation 



(53) 



Similarly 



of' = 



le ) ^-e7r/2 e -i7r/4-(out,int) + e €7r/2 e j7r/4-(out,ext)j ^ 



e-K/2 p iTi/A Wout,ext)^t _|_ p«r/2p-Mr/4^(out,int)^t 



2tt 

6 (m) = ^2 ) J, 



2ir 



a; 



')*]■ 



-(out 



,int) _ nl+K) 



a 



(out,ext) 



2tt 
ra + ie) 



eir/2 p -i7r/4^(in)^t _|_ p -en /2 p in / '4~(in) 



e ' e 



e ' a' 



2tt 



eCT /2 e - l7r /4^(in) + e -«r/2 e i W /4^in))t 



(54) 
(55) 



(56) 
(57) 
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From the Bolgoluibov transformation and the mathematical characterization of |0, in) we 
find that 

|0, in) = iVexp jijT e" |e| 

where AMs a normalization factor. We have therefore exhibited the physical ground state 
of particle-hole pairs correlated between the interior and exterior regions. We now 
have the same formal situation as described in [T2] . If we trace out the "unobservable" 
interior of the black hole, we end up with density matrix is of the form 

p = J2 e ~ 2nlt{i)l \ i > ext ) ® (*,ext|, ( 59 ) 

i 

where % labels the many-body state whose energy is e(i). However, unlike the situation in the 
Unruh-Rindler vacuum [TU [12] our system contains genuine radiation rather that a thermal 
bath. This is because the chiral character of the particles means that they can only flow 
outwards. 



j^(out,ext)\t^(out,int) _|_ Q (out,ext) Wout,int)\f 



at 



2tt 



MO, out), (58) 



IV. DISCUSSION 

The effective space-time metric in which the chiral edge-mode fermions move is 

ds 2 = „ 1 . x dy 2 - dt 2 (60) 
v cd g M 

The quantization of chiral fermions in such a background metric with general v e d ge {y) has 
been carried out in [IB] , although these authors did not consider the effect of an event 
horizon. 

In our case f e d ge = &y, k = X/eB, and a change to an exterior tortoise co-ordinate 



y* = k 1 ln(y) in (60) leads to 

ds 2 = dy 2 - dt 2 . (61) 

The new coordinates reveal that our space-time is flat, but the singularity at the horizon 
is not removed. It has been pushed to y* = — oo, and the interior of the black hole has 
become invisible. A superfluid system with this metric and event horizon was studied by 
Volovik in [TU]. He uses a WKB analytic continuation method to compute the Bogoliubov 
coefficients, and finds the same Hawking temperature as our present calculation, but his 
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non-chiral system has no actual radiation. The agreement in the temperature is perhaps not 
surprising. It must be obvious from looking at the classical trajectories of our particles that 
there is some connection between our 2DEG problem and that of Landau- Zener tunneling 
through an avoided level crossing. Indeed, although the physics is superficially different, 
the Landau-Zener time-dependent Schrodinger equation is solved using the same families of 
parabolic cylinder functions that we have used [17] , and it is well known that an analytically 
continued form of the WKB approximation obtains the correct asymptotic Landau-Zener 
tunnelling probabilities [T8] . 

The most remarkable property of the present model is that the emitted radiation is exactly 
thermal. There is no immediately obvious reason why the mathematical properties of the 
parabolic cylinder functions should lead to this result. In a real black hole the emitted 
radiation is modified by grey-body factors in dimensions greater than two, but that the 
hole can only be in equilibrium with radiation at Xfiawking follows from the geometry of the 
Euclidean section of space-time being asymptotically periodic in imaginary time |19j . Does 
our space-time geometry tacitly force a Euclidean temporal periodicity? 

We can write 

ds* = JL (df _ ^ d{K tf) (62) 
ft y 

so, up to a conformal factor K,~ 2 y~ 2 , the metric is that of Rindler space whose Euclidean 
section t i-> ir has metric 

^Rindler = V^Kt) 2 + dy 2 . (63) 



The absence of a conical singularity at y — in the manifold described by (63) requires 
identifying kt ~ kt + 2ir and so implies a temperature T = 1%k/2tx — which is exactly what 
the tunneling calculation gives. However, given that it blows up at the point of interest, it 
seems unreasonable to ignore the conformal factor, making this argument at most suggestive. 
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